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Abstract, beginabstract: Command not found. 

Adopting ne^w s-process nucleosynthesis scenario and 
branch s-process path, we calculate the heavy-element 
abundances of solar metallicity 3Mq thermal pulse AGB 
(hereafter TP-AGB) stars, and then discuss the correla- 
tion bet"ween heavy-element abundances and C/0 ratio. 
^'^C(a,n)^^0 reaction is the major neutron source, ■which 
is released in radiative condition during the interpulse pe- 
riod, hence gives rise to an efficient s-processing that de- 
pends on the ^^C profile in the ^'^C pocket. A second small 
neutron burst from the ^^Ne source marginally operates 
during convective pulses over previously s-processed mate- 
rial diluted -with fresh Fe seed and H-burning ashes. The 
calculated heavy-element abundances and C/0 ratio on 
the surfaces of AGB stars are compared -with the obser- 
vations of MS, S and C (N-type) stars. The observations 
are characterized by a spread in neutron exposures: 0.5- 
2.5 times of the corresponding exposures reached in the 
three zones of the ^'^C profile sho^wed by Fig. 1 of Gallino 
et al. (1998). The evolutionary sequence from M to S to 
C stars is explained naturally by the calculated heavy- 
element abundances and C/0 ratio. 

Then the heavy-element abundances on the surfaces 
of TP-AGB stars are used to calculate the heavy-element 
overabundances of barium stars, ■which are generally be- 
lieved to belong to binary systems and their heavy-element 
overabundances are produced by the accreting material 
from the companions (the former TP-AGB stars and the 
present ■white d^warfs). To achieve this, firstly, the change 
equations of binary orbital elements are obtained by tak- 
ing the angular momentum conservation in place of the 
tangential momentum conservation, and considering the 
change of Sr/r term; then the heavy-element overabun- 
dances of barium stars are calculated, in a self-consistent 
manner, through ■wind accretion during successive pulsed 
mass ejection, followed by mixing. The calculated rela- 
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tionships of heavy-element abundances to orbital periods 
P of barium stars can fit the observations within the error 
ranges. Moreover, the calculated abundances of nuclei of 
different atomic charge Z, corresponding to different neu- 
tron exposures of TP-AGB stars, can fit the observational 
heavy-element abundances of 14 barium stars in the er- 
ror ranges. Our results suggest that the barium stars with 
longer orbital period P > 1600 d may form through ac- 
creting part of the ejecta from the intrinsic AGB stars 
through stellar wind, and the mass accretion rate is in 
the range of 0.1-0.5 times of Bondi-Hoyle's accretion rate. 
Those with shorter orbital period P < 600 d may be 
formed through dynamically stable late case C mass trans- 
fer or common envelope ejection. 

Key -words: Nuclear reactions, nucleosynthesis, 
abundances- Stars: Abundances- Stars: AGB- Stars: 
Carbon- Stars: Mass loss- Stars: Binaries 



1. Introduction 

Since Burbidge et al. (1957) and Cameron (1957) pub- 
lished their creative studies, the nucleosynthesis theory 
has been developed in deep degree. In particular, asymp- 
totic giant branch (hereafter AGB) stars are very im- 
portant to study element nucleosynthesis and the Galac- 
tic chemical evolution because they synthesize significant 
parts of slow process (hereafter s-process) neutron capture 
elements and ^^C. The products are taken out from the 
stellar interior, He-intershell, to the surface by the third 
dredge-up (hereafter TDU) process, and then are ejected 
into interstellar medium with the progressive stellar wind 
mass loss. 

Our understanding of the AGB nucleosynthesis has un- 
dergone major revisions in these years. The earlier studies 
(Iben 1975; Truran & Iben 1977) illustrated that inter- 
mediate mass TP-AGB stars with ^^Ne(a, n)^^Mg neu- 
tron source (the typical neutron density is least of a or- 



2 



Y. C. Liang et al.: heavy-element abundances of AGB stars and formation of Ba stars 



der of 10^ — 10^°n cm^'^) arc the suitable nucleosynthe- 
sis sites of s-process elements. But new observations shed 
doubts on the above idea (Busso et al. 1995 and references 
therein; Busso ct al. 1999 and references therein). Busso 
et al. (1995) and Lambert et al. (1995) demonstrated that 
the measured abundances of Rb/Sr, the products of the 
branch in the s-proccss path at ^'''Kr, imply a definitely 
lower neutron density (typical of the order of lO^n cm~^), 
which can be provided via the reaction ^^C(a,n)^^0 at 
low temperature in the He-intershell of low mass AGB 
stars. 

Iben & Renzini (1982a,b) indicated that a suitable 

mechanism operated in low mass stars of low mctallic- 
ity to allow the formation of a semiconvective layer, hence 
the ^^C pocket. The pocket is engulfed by the next convec- 
tive pulse where '^"^C nuclei easily capture a nuclei, release 
neutrons. HoUowell & Iben (1988) confirmed the possi- 
bility of formation of a consistent ^^C pocket through a 
local time-dependent treatment of semiconvection. How- 
ever, the semiconvection mixing mechanism was not found 
to work for the ^^C-enriched Population I red giants, like 
the peculiar stars of MS, S and C stars, showing overabun- 
dances of s-process elements in their spectra. 

Straniero et al. (1995) investigated the effect of a pos- 
sible mixing of protons into a thin zone at the top of the 
carbon-rich region during each dredge-up episode, hence 
the formation of ^^C pocket. They suggested that the ^^C 
was completely burnt in the radiative condition, and the 
resulting s-process nucleosynthesis occurs during the qui- 
escent interpulse period, instead of the convective thermal 
pulse. ^^Nc(q!, n)^'''Mg was still active for a very short pe- 
riod during the convective pulse with minor influence on 
the whole nucleosynthesis. Herwig et al. (1997) and Her- 
wig et al. (1998) supported the formation of ^^C pocket 
via hydrodynamical calculations. 

Straniero et al. (1997) adopted the above new s-process 
nucleosynthesis scenario to calculate the s-process nucle- 
osynthesis of solar metallicity low mass AGB stars with 
1 < M/Mq < 3, and gave the detailed results. 

Recently, Gallino et al. (1998) explained further and 
developed the aforesaid new scenario. They divided the 
^^C pocket, q layer, into three zones in the light of the dis- 
tribution in the mass of hydrogen introduced in the ^^C- 
rich intershell. The characteristic neutron exposures in the 
three layers are different. Moreover, when the nucleosyn- 
thesis occurs in a radiative layer, only the nucleosynthesis 
products are ingested into the convective thermal pulse, 
which makes the classical concept of mean neutron expo- 
sure (tq) become meaningless and the simple assumption 
of an exponential distribution of the neutron exposure fail 
to account for the complexity of the phenomenon (Arlan- 
dini et al. 1995, Gallino et al. 1998). Busso et al. (1999) 
reviewed this new s-process nucleosynthesis scenario in de- 
tails. 

The spectral and luminosity studies of AGB stars (in- 
cluding MS, S and N-type C stars) have shown that 



the M^S^C sequence is the result that the low-mass 
AGB stars have undergone carbon synthesis, s-process nu- 
cleosynthesis and the third dredge-ups (Lambert 1991). 
These stars are in the course of experiencing thermal pulse 
stage, and the original chemical abundances of their atmo- 
spheric envelope have been modified by two mixing mecha- 
nisms, namely, the first dredge-up when they became red 
giants and the third dredge-up when they became TP- 
AGB stars (Boothroyd & Sackmann 1988a, b, c and d). 

The predicted evolutionary sequence of M^S^C in 
the heavier— lighter s-element abundances relationship 
(here and hereafter, 'heavier' refers to the second metal 
peak elements: Ba, La, Ce, Nd and Sm etc.; 'lighter' 
refers to the first metal peak elements: Y, Zr etc.) and 
the heavy-element abundances— C/0 relationship (^^C, 
together with the s-process elements, is dredged up from 
stellar interior during the third dredge-up) are important 
to understand the nucleosynthesis and evolution of AGB 
stars. Because (1) they will provide theoretical basis for 
the observed evolution of the sequence, (2) they can check 
the available theories on the evolutions of AGB stars (e.g., 
the beginning of the third dredge- up process, the mass 
and the chemical abundance of the dredge-up material, 
the theory of s-process nucleosynthesis, the stellar wind 
mass loss, and the formation of carbon stars etc.). 

Busso et al. (1992) discussed the heavy-element abun- 
dances of M, MS and S stars using the thermal pulse AGB 
model. Busso et al. (1995) analyzed the heavy-element 
overabundances of carbon stars under the assumption that 
the dredge-up started after reaching the asymptotic dis- 
tribution (about the 20th pulse) . It is difficult to calculate 
the AGB stars evolution and s-process nucleosynthesis. So 
there are few theoretical results to explain the M^S— >C 
evolutionary sequence based on the combination of the 
heavier— lighter s-element abundances ratio and the C/0 
ratio, though the observational abundances exhibit a cer- 
tain regularity. Zhang et al. (1998a) calculated the evolu- 
tion of the surface heavy element abundances and C/0 ra- 
tio for a 3Mq TP-AGB star with initial metallicity 0.015, 
and gave interesting results. But they adopted mean neu- 
tron exposure tq and unbranch s-process path, which have 
been revised in these years. 

In the first part of this paper, we adopt the new s- 
process nucleosynthesis scenario (Straniero et al. 1995; 
Straniero et al. 1997; GalHno et al. 1998; Busso et al. 1999 
etc.), and the branch s-process nucleosynthesis path to 
calculate the s-process nucleosynthesis of solar metallicity 
3Mq AGB stars. And then, we discuss the M^S->C se- 
quence based on the heavy-element abundances and C/0 
ratio. 

The importance of AGB stars nucleosynthesis is not 

only to explain the observational M^S^C sequence but 
also to be responsible for the origin of some other classes 
of stars with overabundances of heavy-elements. 

Observations revealed that some stars with overabun- 
dances of heavy-elements were not luminous to up to the 
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stage of AGB. Following Lambert (1991), the stars show- 
ing heavy-element overabundances are divided into two 
classes: (1) intrinsic TP-AGB stars— they include MS, S 
and C (N-typc) stars exhibiting the unstable nucleus ^^Tc 
(n = 2 X 10° yr) as evidence that they are presently 
undergoing nucleosynthesis activity and the third dredge- 
up, and (2) extrinsic AGB stars— they include the various 
classes of G-, K-type barium stars and the cooler S, C 
stars where ^^Tc is not observed. It is generally believed 
that the extrinsic AGB stars belong to binary systems and 
their heavy-element overabundances come from accretion 
of the matter ejected by the companions (the former AGB 
stars, now evolved into white dwarfs) (McClure et al. 1980; 
BofRn & Jorissen 1988; Jorissen et al. 1998; Jorissen & 
Van Eck 2000; Jorissen 1999). The mass exchange took 
place about 1 x 10^ years ago, so the ^^Tc produced in the 
original TP-AGB stars have decayed. The accretion may 
either be disk accretion (Iben & Tutukov 1985) or com- 
mon envelope ejection (Paczynski 1976). Han et al. (1995) 
detailedly investigated the evolutionary channels for the 
formation of barium stars. In this paper we will only dis- 
cuss the stellar wind accretion scenario because it is very 
important to explain the formation of barium stars (Boffin 
& Jorissen 1988; Jorissen et al. 1998). 

BofRn & Jorissen (1988) calculated qualitatively the 
variation of orbital elements caused by wind accretion in 
binary systems. They also estimated the heavy-element 
overabundances of barium stars. Subsequent papers (Zacs 
1994, Boffin & Zacs 1994) used similar methods to calcu- 
late the overabundances, and interpreted the relationship 
between the heavy-element abundances and the orbital 
periods of barium stars. 

Some important conclusions have been drawn in the 
theory of wind accretion, but the previous calculations on 
orbital elements were not very reliable because of the ne- 
glect of the 5r jr term (r represents the distance between 
the two components of the binary system ), and using the 
tangential momentum conservation (Chang et al. 1997 and 
references therein). For the rotating binary system with 
wind mass loss, the total angular momentum conserva- 
tion is more reasonable than tangential momentum con- 
servation. Also, the previous calculations of heavy-element 
overabundances used the 'step-process' (Boffin & Jorissen 
1988, Boffin & Zacs 1994), which means that the over- 
abundance factor changes at one single instant from 1 to 
/ (/ is the relative ratio of the heavy-elements to the solar 
abundances, and differs for different elements. Earlier cal- 
culations used the mean / value of carbon stars) , and then 
keeps the value until the end of the AGB phase. However, 
after the start of the TDU, the overabundance factor / of 
the intrinsic AGB stars changes during successive dredge- 
ups. It is after a number of dredge-ups that the C/0 ratio 
in the outer envelope of the intrinsic AGB star reaches the 
value 1, which means that the AGB star becomes a carbon 
star. The heavy-element overabundances of the barium 



star should be caused by the successive pulsed accretions 
and mixing. 

According to the analysis to the orbital elements of 
barium and S stars, Jorissen & Mayor (1992) presented 
the evolutionary pathways of binaries leading to barium 
and S systems. They concluded that the binary systems 
with longer orbital period formed through wind accretion 
and those with shorter orbital period formed via Roche 
lobe overflow. But the specific range of orbital period was 
not presented. Jorissen et al. (1998) analyzed the orbital 
elements of a large sample of binary systems to give insight 
into the formations of barium and Tc-poor S stars. They 
suggested that barium stars with orbital period P>1500 
days formed through wind accretion scenario. Zhang et al. 
(1999) suggested that the barium stars with P>1600 days 
formed via wind accretion according to their model. Liu 
et al. (2000) discuss the evolutions of orbital elements of 
barium stars from normal G, K giants. 

Besides the orbital elements, the heavy-element abun- 
dances of barium stars have been discussed in some lit- 
eratures. Busso et al. (1995) discussed the observational 
heavy-element abimdances of barium stars. And a more 
detailed analysis of the abundance distributions for five 
stars has been performed using the method of mixing the 
accretion mass with the envelope mass. But the effect of 
mass accretion and the changes of orbital elements were 
not considered. Chang et al. (1997) and Liang et al. (1999) 
attempted to explore the relationships between the heavy- 
element overabundances and orbital elements of bariums 
stars using the binary accretion scenario, but with the 
shortcomings: or using the tangential momentum conser- 
vation, or adopting the old nucleosynthesis scenario of TP- 
AGB stars. 

In the second part of this paper, we firstly reduced the 
variation equations of orbital elements based on the an- 
gular momentum conservation model of wind accretion, 
then we calculate the heavy-element overabimdances of 
barium stars via successive pulsed accreting matter en- 
riched heavy-elements from the intrinsic AGB stars, and 
mixing the matter with their envelopes. 

This paper is organized as follows. The observational 
data of MS, S, C (N-type) stars and barium stars are 
given in Sect. 2. In Sect. 3, we present the model and the 
main parameters of AGB stars nucleosynthesis and the 
angular momentum conservation model of wind accretion 
scenario for barium stars. Sect. 4 illustrates and analyzes 
our results. We conclude and discuss in Sect. 5. 

2. Observations 

The abundance data of the heavy-elements, C and O 

elements of MS and S stars are taken from Smith & 
Lambert (1985, 1986, 1990). The listed values of [N/Fe] 
(=log(N/Fe)star— log(N/Fe)sun) have typical errors 0.1-0.2 
dex, and N represents a kind of s-process element. Indeed, 
for the observations in the series papers of Smith & Lam- 
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bert, the accuracy of the iron peak elements (Ti, Fc and 
Ni) abundances and the heavy-elements Y, Zr, Nd abun- 
dances were better. But the abundance of Nd element is 
only at 50% of s-origin, it is not reasonable to regard Nd 
as the representative of heavier s-elements (Busso et al. 
1992; Zhang et al. 1998a), and the reasonable representa- 
tive should be the mean abundance of Ba, La, Nd and Sm 
(Busso et al. 1995). In order to do this we suggest that the 
initial abundance of Nd have to be increased by 0.11 dex 
before the data of MS and S giants are compared with car- 
bon stars (Busso et al. 1995). Y, Zr are the representatives 
of lighter s-elements. For C and O elements, the accuracy 
is higher at 0.1 dex in Smith & Lambert (1990). Consider- 
ing that, at present, the sample with both observed C/0 
ratios and heavy-element abundances is rather limited and 
there is a close relationship between the Tc-no MS, S stars 
(stars that contain no Tc) and the thermal pulse AGB 
stars (Busso et al. 1992; Jorissen & Mayor 1992), we can 
usefully include the Tc-no stars in our study of evolu- 
tionary tendency of the M^S^C sequence (Busso et al. 
1992). 

The abundance data of heavy-elements of the 12 car- 
bon stars (N-type) are taken from Utsumi (1985), with 
typical errors of 0.4 dex in [N/Ti], which are dealt like in 
Busso ct al. (1995). The corresponding abundances of C, 
O elements are taken from Lambert et al. (1986). 

The heavy-element abundances of 12 barium stars 
arc taken from Zacs (1994). For ( Cap and HR 774, 
the heavy-element abundances are from the same 
sources as in Busso et al. (1995). The sun, instead of 
e Vir, is chosen to be a comparison standard, namely, 
[N/Fe]=log(N/Fe)star-log(N/Fe)sun- The abundances of e 
Vir to the sun are taken from Tomkin & Lambert (1986). 
The new solar system abundances are taken from Anders 
& Grevesse (1989). The overabundances of s-process ele- 
ments in barium stars are characterized by a mean value of 
[s/Fe] = ([Y/Fe] + [Ba/Fe] + [La/Fc] + [Cc/Fc] + [Nd/Fe])/5.0. 
The typical error is 0.2 dex (BofHn & Zacs 1994). Zr has 
been excluded from this mean value, since the abundance 
of Zr is poorly related to other heavy-elements in the 
sample (Boffin & Zacs 1994). The corresponding obser- 
vational data of orbital elements, orbital period P and 
eccentricity e, are taken from McClure & Woodsworth 
(1990), Boffin & Zacs (1994) and Jorissen et al. (1998). 

3. The model and the main parameters 

3.1 The nucleosynthesis and mixing in AGB stars 

AGB star consists of three layers: a C-0 core of degen- 
erate electron, a He and H double burning shell, and a 
convcctivc envelope. Heavy-elements synthesized in the 
He-intershell are dredged up to the convective envelope 
by the convection that appears in the cooling contrac- 
tion after a thermal pulse, and show up as the observed 
overabundances of the heavy-elements and ^^C. With the 



evolution of the star, its core mass steadily increases with 
increasing thermal pulses, while the mass of the convective 
envelope decreases correspondingly via wind mass loss, un- 
til it is completely exhausted, while marks the life span of 
the thermal pulse AGB star. 

During every thermal pulse, we calculate the s-process 
nucleosynthesis. ^^Fe is taken to be the seed nucleus for the 
s-process. The reaction 209Bi(n,7)2i0Bi(a)20'3pb(2i0Bi has 
an a-decay half-life of 5.01 days) terminates the reaction 
chain of the slow process of neutron capture. The branch 
s-process path is adopted in the calculation. The neutron 
capture cross sections of most of nuclei, are taken from 
Beer et al. (1992). The /3-dccay rates and electron capture 
rates of nuclei are taken from Takahashi & Yokoi (1987). 
More recent values of some nuclei are taken from the same 
literatures as those cited by Gallino et al. (1998). For the 
initial abundances, we use standard red giant abundances, 
which differ from the solar abundances only by a constant 
factor: Ni = NiQ x Z/Zq. At present, the solar system 
abundances are the most detailed abundance distribution 
obtained by us. According to the results of Straniero et 
al. (1997) and Gallino et al. (1998), the main parameters 
adopted in calculation are as following. 

(1) Core mass At the onset of the thermal pulse, the 
mass of the C-0 core is 0.572Mq. With successive pulses, 
both the H burning and He burning shells move outwards 
(in mass coordinates) and so does the C-0 core. The core 
mass is up to 0.611Mq at the 8th pulse, during which the 
third dredge- up begins. The variation of the core mass 
influence directly the wind mass loss. The core masses of 
stars in different TP periods are taken from table 4 of 
Straniero et al. (1997). 

(2) Mass loss through wind Wind mass loss is one 
of the most important ingredients in the computation of 
both the nucleosynthesis occurring in TP- AGB stars and 
the wind accretion on the secondary component of a bi- 
nary system. Here, wc adopt the mass loss rates given by 
Straniero et al. (1997) using Reimers formula (Reimers 
1975) based on the TP- AGB stellar models (their table 
4). 

The mass lost during the time Ai, the interpulse pe- 
riod, is 

AM = -MAt, (1) 

which is taken from table 4 of Straniero et al. (1997). 

(3) Dilution Factor The physical process that in- 
fluences signiflcantly the surface abundance is the third 
dredge-up. The dilution factor / is deflned as the ratio of 
the mass dredged in every pulse to the mass of the enve- 
lope after the third dredge-up begins, namely, 

^ AMtdu 
^ Me„v + AMtdu' ^ ' 

where AMtdu is the mass dredged up from He-intershell 

at every pulse since TDU begins, and M^nv is the mass of 
the envelope. The cumulative effect of the various TDU, 
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^AMtdu, have been given in Fig. 9a of Gallino et al. 
(1998). 

(4) Overlap factor The overlap factor rr, rr = 

1 — AMc/McsH, is the ratio of the number of a kind of 
nuclide that is exposed to neutron in two successive pulses 
to the number exposed in the previous pulse, that is, the 
overlap fraction of the convective shells between the two 
successive pulses. Where AM^ is the increment of the core 
mass during the previous interpulse period, and Mcsh is 
the mass of the convective inter-shell. The variation of rr 
with core mass was given by Iben (1977) and Renzini & 
Voli (1981). According to suggestion of Straniero et al. 
(1995), rr steadily decreases with increasing core mass. In 
this paper, we adopt the values given by Fig. 9c of Gallino 
et al. (1998). 

(5) Neutron source -'^■'C(a, n)-'^^0 reaction is the ma- 
jor neutron source of s-process nucleosynthesis, which is 
released in radiative conditions at Tg ~ 0.9 during the in- 
terpulse period, hence gives rise to an efficient s-processing 
that depends on the ^'^C profile in the ^"^C pocket, q layer. 
The layer is divided into three zones, in which the neu- 
tron exposures change with the various interpulse periods 
(Gallino et al. 1998, see their Fig. 6). Then the products 
of nucleosynthesis are ingested into the convective thermal 
pulse, and mixed with s-process material already from the 
previous pulses, and with the H-burning ashes from the 
below H shell. A second small neutron burst from the ^^Ne 
source operates during convective pulses due to the high 
neutron density and high temperature. The ^^Ne source to 
the total neutron irradiation is small, ranging from ^0.002 
mb"-'^ at the 9th pulse to ^^0.05 mb"-'^ in the latest pulses, 
where the average effective temperature is ~23 keV. 

(6) Light neutron poisons Some isotopes can consume 
neutron produced by the neutron source isotopes, ^"'C or 
^^Ne, through some reactions. They are named as neutron 
poisons. For example, the isotope ^^N through its reso- 
nant channel ^''N(n,p)^^C, ^^Al in its long-lived groimd 
state synthesized by the reaction ■^^Mg{p,j)^^A\{p,'^)^'^Si 
(Gallino et al. 1998 and references therein). However, in 
the previous convective pulses, ^''N is destroyed in the 
reaction chain leading to ^^Ne production, and a large 
fraction of ^^Al has undergone substantial depletion by 
neutron capture. The above discussion makes it clear that 
compared to a convective burning scenario, the radiative 
s-process is much less affected by the filtering effect of light 
neutron poisons (Gallino et al. 1998). So in the radiative 
s-process nucleosynthesis calculation, we neglect the influ- 
ence of these light neutron poisons (the slightly effects on 
nucleosynthesis calculation caused by this will be further 
discussed in section 4.1). 

(7) Value of C/0 At the beginning of the AGB 
phase, the C/0 on the surface of a low mass stars is ap- 
proximately 0.31, and lower than the initial (solar) value 
(approximately 0.45) because it has been modified by the 
first dredge- up process, which reduces the ^^C abundance 
of the surface by about 30%. With the third dredge- up 



process, ^^C of He-intershell is mixed into the stellar sur- 
face with the mixing of s-process material. So the C/0 
ratio on the surface increases gradually. In calculation, we 
adopt the C/0 ratio from table 4 of Straniero et al. (1997). 

.3.2 The heavy- element overabundances of barium stars 

3.2.1 The angular momentum conservation model of wind 
accretion 

For the binary system, the two components (an intrinsic 
AGB star, the present white dwarf, with mass Mi, and 
a main sequence star, the present barium star, with mass 
M2) rotating around the mass core C, so the total angu- 
lar momentum is conservative in the mass core reference 
frame. If the two components exchange material through 
wind accretion, the angular momentum conservation of 
total system is showed by: 

A{fj,r'^e) = wrf (AMi + AM2) + r2v{AMi + AM2), (3) 

where /x is reduced mass, and r is the distance from M2 to 

Ml . ri, r2 are the distances from Mi , M2 to the mass core 
C respectively, w (= 2n/P, where P=2TrA'^{l - e^)^ /h 
is orbital period) is angular velocity, v is an additional 
effective velocity defined through the angular momentum 
variation in the direction of orbital motion of component 
2. The first term on the right side of the equal-sign is the 
angular momentum lost by the escaping material, and the 
second term is the additional angular momentum lost by 
the escaping material. 

Using the similar method to that adopted by Huang 
(1956), Boffin & Jorissen (1988) and Theuns et al. (1996), 
considering the angular momentum conservation of total 
system and not neglecting the square and higher power 
terms of eccentricity, we can obtain the change equations 
of the orbital elements: 



AA 



= -2(1 -e^) 
+ 2(1 -e^) 
+ [2(l-e2)5 



AMi AM2 AMi + AM2 



Ml 
M2{AMi 



M2 
AM2) 



Mo 



Vorh 



Mi{Mi+M2) 

. AMi + AM2 



1]- 



M1 + M2 



(4) 



eAe 



= [l-(l-e^)^ 



AMi AM2 AMi + AM2 



Ml 



M2 



Ml + M2 



- [l-(l-e^) 



2 i^M2(AMi -I- AM2) 



Mi(Mi+M2) 
3e2 AMi + AM2 V 



2(l-e2)5 



Ma 



forb 



(5) 



where A is the semi-major axis of the relative orbit 
of component 2 around 1, and e is the eccentricity (more 
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details can be found in Appendix and Liu et al. (2000)). 
Here, we take v = (BofHn & Jorissen 1988). 

For the mass accreted by the barium star, we use 
the Bondi-Hoyle (hereafter B-H) accretion rate (Bondi & 
Hoyle 1944; Theuns et al. 1996; Jorissen et al. 1998): 



1 + ('J^orb/'f^ej)^ 



AMi, (6) 



where a is a constant expressing the accretion efficiency, 
and its value is taken from Theuns et al. (1996). Theuns 
et al. (1996) indicated that taking a between 0.5 and 1, 
as suggested by the numerical simulations of Ruffert & 
Arnctt (1994), the actual accretion rate deduced from the 
smoothed particle hydrodynamics (hereafter SPH) simu- 
lation in the 7 = 1.5 case was thus about 10 times smaller 
than that predicted by the B-H formula. Here we take 
a=l for the B-H formula. In fact. Boffin & Zacs (1994) 
suggested that the actual accretion rate is between 0.1 
and 1 times of B-H rate. We take 0.15 times of the B-H 
rate as the actual accretion rate for our standard case. Vej 
is the wind velocity, and Voib is the orbital mean veloc- 
ity. After fixing the initial conditions, for the mass AMi, 
ejected at each pulse by the primary star, we can solve 
the Eqs.(4)-(6) for AM2, the accreted mass by the sec- 
ondary star. And then, the heavy-element abundances on 
the surface of barium stars can be calculated. 



3.2.2 The heavy-element overabundances of barium stars 

The calculation is completed by two separated steps. 
Firstly, adopting the theory of s-process nucleosynthesis 
and the latest TP- AGB model, we calculate the degree of 
the overabundances of the intrinsic AGB star (the present 
white dwarf) at each ejection. Then, combining the accret- 
ing matter predicted by the model of wind accretion on 
successive occasions and mixing, we calculate the heavy- 
element overabundances of the barium star. The over- 
abundance factor for nuclide i on the secondary star, Qi, 
is given by 



9i = 



n=l,m 
n=l.m 



(7) 



where M^nv is the mass of the outer envelope of barium 
star, M^cc is the mass accreted by barium star during the 
period of the n-th ejection of the intrinsic star, and /" is 
the overabundance factor of the nuclide i of the intrinsic 
AGB star on that occasion, m is the total ejecting number 
undergone by the intrinsic AGB star, is the overabun- 
dance factor of nuclide i of barium star before the mass 
accretion. In the above formula we have assumed that the 
accreting matter has already completely mixed with the 
outer convective envelope of barium star. 
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Fig. 1. Comparison of the theoretical predictions of sur- 
face heavy-element abundances for different neutron ex- 
posures (solid lines) with observations of intrinsic AGB 
stars. The numbers of 'a' represent the times of the cor- 
responding exposures in the ^^C profile suggested by the 
Fig. 1 of GaUino et al. (1998). 



We take as standard case: Mi,o=3.OM0, M2,o=l-3M0, 
t;ej=15 km s~^ and 0.15 times of the Bondi-Hoyle's accre- 
tion rate. 

4. Results and analysis 

4-1 The intrinsic AGB stars 

Using the recent evolutionary model and nucleosynthesis 
scenario of intrinsic TP-AGB stars (Straniero et al. 1995; 
Straniero et al. 1997; GaUino et al. 1998), we calculate the 
heavy-element abundances on the surface of solar metallic- 
ity 3Mq AGB star with wind mass loss. The TDU begins 
at the 8th pulse. The average heavy element abundances 
on the surface are given in Fig. 1. 

In Fig. 1, the abscissa represents the average logarith- 
mic ratio of the abimdanccs of the heavier s-clcmcnts Ba, 
La, Nd, and Sm ('hs') to the lighter s-elements Y and Zr 
('Is'), which depends on the neutron exposure. The ordi- 
nate is the logarithmic enhancement in 'Is' with respect 
to iron in the envelope, which is closely related to the 
synthesis in the He-intershell and the dilution factor. The 
different curves represent the results of different neutron 
exposures. 

Indeed, the neutron exposures of ^^Ne in different ther- 
mal pulses have been showed in Sect. 3.1. In calculation, 
we adjust the neutron exposure of ^^C source. GaUino et 
al (1998) suggested the ^^C profile in ^^C pocket by their 
Fig. 1, namely their 'standard' case. Then they showed 
that the resulting heavy-elements abundances were non- 
solar by their Fig. 14, while simply increasing the previ- 
ously adopted abundance of -'^^C by a factor of 2 could 
reproduce the main component of solar system (see their 
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Fig. 16). Here, instead of the abundance of ""^"^C, we ad- 
just the total neutron exposure caused by ^'^C source to 
calculate the nucleosynthesis of AGB stars. Thus, the dif- 
ferent values of 'a' in our Fig. 1 represent the different 
times of the actual total neutron exposures to the 'stan- 
dard' case suggested by Gallino et al. (1998). Our result of 
a=1.0 corresponds to the 'standard' case of Gallino et al. 
(1998), and the a=1.5 case corresponds to their results of 
increasing the previously adopted abundance of ^^C by a 
factor of 2. In reality, Gallino et al. (1998) have illustrated 
this times relationship (1.5 vs. 2) in the description and 
discussion about their Figs. 14, 15 and 16. 

Our Fig. 1 shows that our results arc very similar to 
those of Gallino et al. (1998). The curve with a=1.0 ex- 
hibits clearly that the abundances of the heavier nuclei 
arc lower than those of the lighter nuclei in the 23rd inter- 
pulse phase like Galhno et al. (1998), namely [hs/ls]<0.0 
(the abscissa), which means that this case can not repro- 
duce the solar abundance distribution due to a lower pro- 
duction for heavier nuclei. While the curve with a=1.5 
produces [hs/ls]~0.0 in the 23rd interpulse phase, which 
means that the abundances of the heavier and the lighter 
nuclei are solar like. So, our 'average' heavy-element abun- 
dances on the surfaces of AGB stars can fit the detailed 
abundance distribution of nuclei obtained by Gallino et 
al. (1998). In addition, we should note that, because we 
neglect the effects of light neutron poisons on nucleosyn- 
thesis in calculation, our 'a' values are slightly higher than 
the corresponding case of Gallino et al. (1998). 

We calculate some other curves with different a values, 
which arc given in Fig. 1 too. The results show that all of 
our available observational data are compatible with the 
range 0.5<a<2.0. The evolutionary curves for three cases 
(a=1.0, 1.5, 2.0) move upward as the drcdge-ups proceed, 
and from the region of MS, S stars to reach the region of 
the carbon stars (symbol plus). And the larger neutron 
exposure, the higher [hs/ls] ratio will be, which means 
the large neutron exposure benefits the production of the 
heavier s-process elements. 

^^C isotope in the Hc-intershell of TP- AGB stars, to- 
gether with the s-process elements, is dredged-up and 
mixed to the envelope by the TDU process, which causes 
that the C/0 ratio on the surface increases gradually with 
the thermal pulses. The increase correlates to the over- 
abundances of heavy-elements. Our Fig. 2 displays the 
relationship of heavy-element abundances to C/0 ratio. 
The values of 'a' are relevant to the values in Fig. 1. 
From the beginning of TDU, the values of C/0 ratio and 
heavy-element abundances increase gradually. After some 
dredge-ups, C/0 becomes greater than 1, that is, the stars 
become carbon stars (a=1.2, 1.5, 2.5). The different curves 
show that the higher neutron exposure, the more heavy el- 
ements are produced. 

The results of Figs. 1 and 2 illustrate that, adopting the 
new TP- AGB nucleosynthesis scenario, choosing the rea- 
sonable parameters, we can explain the observed heavy- 



Z J 


. + C (N-type) 










^ MS,S (Tc-yes) 






2.5 




^ MS,S (Tc-no) 






A ' 








1-5 


Li_ 








1.2 
1.0 




















0.8 


























^-^^'^^ 




0.5 _ 


0.0 











0.0 0.5 1.0 1.5 

C/0 



Fig. 2. Comparison of theoretical predictions of the re- 
lationships between heavy-element abundances and C/0 
ratios for different neutron exposures (solid lines) with 
observations of intrinsic AGB stars. The numbers of 
'a'represent the same meanings as in Fig. 1. 

element abundances and C/0 ratio of MS, S and C (N- 
type) stars. Also, we can explain the M— >S— >C evolu- 
tionary sequence on the basis of the lighter— heavier s- 
element abundances relationship and the heavy-element 
abundances— C/0 ratio relationship simultaneously. 

^.2 The barium stars 

Our model assumes that the binary system remains al- 
ways detached, so that the only interaction between the 
two components of the binary system occurs via the ac- 
cretion by the less evolved component of some fraction 
of the mass lost through stellar wind by the carbon- and 
heavy elements-rich AGB component. In other words, the 
internal structure of each component remains unaffected 
by the presence of its companion, so that the usual struc- 
tural properties of single stars remain valid. Based the 
above-mentioned model, the heavy-element abundances 
of intrinsic AGB stars are calculated firstly, then the 
heavy-element overabundances of barium stars are self- 
consistently calculated through the progressive pulsed 
wind accretion and mixing. 

The change equations of orbital elements are ob- 
tained by adopting the angular momentum conservation 
of the total system (including the two companions and the 
ejected matter) (see Sect. 3.2.1 and Appendix). 

The observed orbital periods of 14 barium stars are in 
the range of 80.53 to 6489 days. According to the discus- 
sions of Jorissen et al. (1998) and Zhang et al. (1999), the 
barium stars with orbital period P > 1500 or P > 1600 
days formed through wind accretion, while those with 
P < 600 days formed through other scenarios. So in calcu- 
lation, we adopt 600-7000 d to be the orbital period range 
of wind accretion. 
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Fig. 3. Comparison of the predicted to observed relation- 
ships of surface heavy-element abundances to the orbital 
periods of barium stars corresponding to different neutron 
exposures of intrinsic AGB stars (solid lines) in standard 
case of wind accretion. The numbers represent the same 
meanings as in Fig. 1. 

The actual mass accretion rate between the compo- 
nents of binary systems can be 0.1-1 times of B-H rate 
(Boffin & Zacs 1994), while SPH simulation (Theuns et al. 
1996) indicated that the actual rate was about 10 times 
smaller than B-H rate. We adopt 0.15 times of B-H rate 
as the standard case in our calculation. 

The calculated relationships between the heavy- 
clement overabundances [s/Fc] and orbital period P of 
barium stars are displayed in Fig. 3, which are based on 
the standard case of wind accretion. Here, again, the var- 
ious curves correspond to different times of the neutron 
exposures of the intrinsic AGB stars, and 'a' represents 
the same meanings as in Fig. 1 and Fig. 2. Every point 
of every curve refers to the different orbital period, within 
the range 600-7000 d. Since the shorter the orbital pe- 
riod, the greater the accretion, and hence the larger the 
heavy-element overabundances will be. For the curves in 
Fig. 3, the higher points correspond to the shorter peri- 
ods. The theoretical results with 0.8 < a < 2.5 can fit the 
observations within the error range. 

To make further examination to the dependence of 
heavy-element abundances on the angular momentum 
conservation model of wind accretion, the mass accretion 
scenario of barium stars is analyzed carefully by compar- 
ing the predicted heavy-element overabundances of differ- 
ent atomic charge Z with the observations of 14 barium 
stars. 

In calculation, we try to make the calculated eccentric- 
ity e and orbital period P match to the observations of the 
barium stars. The results are given in Figs. 4a — 4n. The 
corresponding neutron exposure of intrinsic AGB compo- 
nents, a, and the orbital period P of barium stars are 
exhibited in the every Figure. 



These results show that, for the 9 long-period barium 
stars (P>1600 d), the calculated curves can fit well to the 
observed heavy-element abundances in the error range (see 
Figs. 4a — 4i) according to the standard case of our wind 
accretion scenario. For the two classical barium giant stars 
HD 204075 (C Cap) and HD 16458 (HR 774), the results 
will fit the observations better after the mass accretion 
rate is improved to 0.5 times as much as the B-H accre- 
tion rate (see Figs.4j, 4fc). For the three barium stars with 
shorter orbital periods, HD 199939, HD 46407 and HD 
77247 (the orbital periods are 584.9, 457.4 and 80.53 days 
respectively), the calculations can not fit the observations 
(sec Figs. 4:1 — n). 

Because the same program barium stars are chosen in 
Fig. 4 as in Fig. 3, the reasonable result would be that 
the same observational sample should correspond to the 
same neutron exposure of intrinsic AGB stars in the two 
figures. The results are advantageous to the standpoint: 
for most of barium stars, the corresponding neutron ex- 
posures of 'a' in Figs. 4a — 4n are consistent to the values of 
'a' in Fig. 3 within the error ranges. Namely, the compar- 
isons between the predicted and the observed results are 
basically identical in the heavy-element abundances dia- 
gram and the heavy-element abundance— orbital period P 
diagram of barium stars. 

Comparisons between the predicted and the observed 
abundances of barium stars in Fig. 3 and Fig. 4 show 
that the barium stars with longer orbital period (P>1600 
d) form through wind accretion, and the change range of 
mass accretion rate should be 0.1 to 0.5 times as much as 
the Bondi-Hoyle's accretion rate; the barium stars with 
shorter period (P<600 d) can form through other scenar- 
ios: stable case C disk accretion or common envelope ejec- 
tion. Our results support quantitatively the conclusions of 
Jorissen et al. (1998) and Zhang et al. (1999). 

5. Conclusion and discussion 

5.1 The intrinsic AGB stars 

Adopting the latest AGB stars evolutionary theory 
and nucleosynthesis scenario, the heavy-element abun- 
dances of solar mctallicity 3Mq AGB stars arc cal- 
culated. It is shown that, adopting reasonable param- 
eters, the calculated results of the heavier /lighter s- 
elements ratio— overabundance relationships and the C/0 
ratio— overabundance relationships can fit the observa- 
tions. The evolution of AGB stars along the M^S— >C 
sequence is thus further explained from the evolutionary 
theories and heavy-elements nucleosynthesis of AGB stars. 
A 3Mq AGB star will undergo about 27 pulses before it 
becomes carbon star. 

The results showed that the overabundances of s- 
process elements depend significantly on the neutron ex- 
posure: for the lower neutron exposures (a<1.5), the 
lighter s-process elements are more abundant than the 



9 



heavier ([hs/ls]<0.0); for the higher neutron exposures 
(a>1.5), the latter can increase strongly with the increase 
of neutron exposures, and the lighter s-process elements 
change slightly ([hs/ls]>0.0); while the results with a=1.5 
can explain the solar heavy-element abundance distribu- 
tion, which confirm the results of Gallino et al. (1998). 
Moreover, ^^C abundance correlates to the s-proccss el- 
ement abundances. Both C/0 ratio and s-element abun- 
dances increase with the occurrence of TDU. After some 
dredge-ups, C/0 ratio reaches 1, from which the AGB 
stars become carbon stars. The observed MS, S and C (N- 
type) stars lie in a range of neutron exposure: a=0.5-2.5. 
The over low neutron exposure can not produce carbon 
stars (e.g. a=0.5). Our calculation thus provides a fur- 
ther theoretical basis for the evolution of AGB stars along 
the M^S^C sequence based on the heavy-element abun- 
dances and C/0 ratio. The general agreements between 
our calculations and the observations indicate that the 
parameters, the theories of evohition and nucleosynthesis 
scenario of AGB stars adopted in calculation are reason- 
able. 

The stars with initial main sequence mass 1.5— AMq 
can form carbon stars (Groenewegen et al. 1995). We con- 
sider only a 3Mq star model in our calculation. So the 
comparison between the calculated results and the obser- 
vations is somewhat rough. With further studies of the 
evolutionary theory of AGB stars and increasing observa- 
tional data, we can expect a deeper understanding on the 
nucleosynthesis of AGB stars. 

5.2 The barium stars 

Taking the conservation of angular momentum in place of 

the conservation of tangential momentum for wind accre- 
tion scenario, considering the change of dr/r term, and 
not neglecting the square and higher power terms of ec- 
centricity, the change equations of orbital elements are 
recalculated. We combine wind accretion with the nucle- 
osynthesis of intrinsic AGB stars, to calculate, in a self- 
consistent manner, the heavy-element overabundances of 
barium stars through mass accretion during successive 
pulsed ejection, followed by mixing. 

The calculated relationships of heavy-clement abun- 
dances — orbital period P can fit the observations within 
the error ranges. Moreover, the predictions of the detailed 
abundances of different atomic charge Z can match well 
the observations of 11 program barium stars with longer 
orbital period (P>1600 d). The corresponding neutron 
exposures arc in the range of 0.8<a<2.5. The higher neu- 
tron exposure (e.g. a=2.0) will produce the more abun- 
dant the heavier s-elements than the lighter (see Fig. 4c); 
on the contrary, the abundances of the lighter s-elements 
are higher than the heavier with low neutron exposure 
(e.g. a=1.0, see Fig. 4a); while a=1.5 will produce the 
almost equal abundances of the heavier to the lighter s- 
elements (see Fig. 4h). These neutron exposures can fit 



the corresponding results of intrinsic AGB stars. Natu- 
rally, we can understand the observations of no-Tc MS 
and S stars, which are commonly believed to be the de- 
scendants of barium stars. These results of the extrinsic 
AGB stars confirm the reliability of the nucleosynthesis 
and evolution of intrinsic AGB stars. Simultaneously, the 
results confirm that our wind accretion model and param- 
eters adopted are suitable. 

Analyzing our results, we understand that the barium 
stars with longer orbital period (P>1600 d) form through 
wind accretion. Those with shorter orbital period (P<600 
d) form through other scenarios, such as dynamically sta- 
ble late case C mass transfer or common envelope ejection. 
Moreover, the change range of mass accretion rate should 
be 0.1 to 0.5 times as much as the Bondi-Hoyle's accretion 
rate. The corresponding range of orbital periods and mass 
accretion rate to the formation of barium stars still need 
to be tested by more observations. 

At present, the orbital elements of a large sample bar- 
ium stars, Tc-poor S stars have been published (Udry et 
al. 1998a, 1998b; Carquillat et al. 1998). But the corre- 
sponding heavy-element abundances have not been ob- 
tained. So we need the high resolution, high signal-to-noise 
spectral observations of these stars, which are combined 
with the observations of orbital elements, to research their 
characters and formation. 

In addition, we should note that metallicity is an im- 
portant factor to the AGB stars nucleosynthesis and the 
formation of chemical peculiar stars. Actually, ^^C neu- 
tron source is related to the metallicity (Busso et al. 1995, 
1999 and references therein; Gallino et al. 1999). Gallino 
et al. (1998) calculated the s-element nucleosynthesis of 
2Mq AGB stars with low metallicity Z=0.01, and ob- 
tained similar abundance distribution to the 3Mq with 
solar metallicity AGB stars model. Busso et al. (1999) and 
Zhang et al. (1998b) have discussed the inverse correlation 
between the heavy-element abundances and the metallic- 
ity [Fe/H] of intrinsic and extrinsic AGB stars. Also, the 
nucleosynthesis results of low metallicity AGB stars are 
more suitable to study the Galactic chemical evolution in 
the early stage of the Galaxy. For extrinsic AGB stars, 
Jorissen et al. (1998) suggested that the s-process opera- 
tion was more efficient in a low-metallicity population, so 
the Pop. II CH stars may has accreted the material much 
enriched in heavy elements from the former AGB com- 
panion. In this paper, our main aims are (1) calculating 
the AGB stars nucleosynthesis, so that we can explain the 
observed heavy-element abundances of MS, S and C (N- 
type) stars, which are near solar metallicity, and supply 
evidence to the M^S^C evolutionary sequence; (2) dis- 
cussing the wind accretion scenario of Ba stars, which are 
near solar metallicity too (Zacs 1994; Smith et al. 1993). 
So we only calculate the solar metallicity case. We will ex- 
tend to study the low metallicity case in the forthcoming 
paper. 
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Appendix 

The angular momentum conservation model of 
wind accretion for barium stars: Equations 

For the binary system, the two components (an intrin- 
sic AGB star witli mass Mi, the present white dwarf, and 
a main sequence star with mass M2, the present barium 
star) rotating around the mass core C, so the total angu- 
lar momentum is conservative in the mass core reference 
frame. If the two components exchange material through 
wind accretion, the angular momentum conservation of 
total system is showed by: 

A{ij,r^O) = LorliAMi + AM2) + r2v{AMi + AM2), (8) 

where /z is reduced mass, and r is the distance from M2 to 
Ml. ri and r2 arc the distances from Mi, M2 to the mass 
core C respectively, iv {=2tt/P) is angular velocity, where 
P = 27rA2(l - e'^)i/h is the orbital period (Huang 1956). 
V is an additional effective velocity defined through the 
angular momentum variation in the direction of orbital 
motion of component 2. The first term on the right side 
of the equal-sign is the angular momentum lost by the 
escaping material and the second term is the additional 
angular momentum lost by the escaping material. 

For the binary system, according to Huang (1956), the 
changes of orbital elements, the orbital semi-major axis A 
and eccentricity e, are 
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According to the angular momentum conservation 
model, the A{rO) term can be obtained from the following 
equation: 



M1M2 
Ml -I-M2 



rA{re) = A 



A 



M1M2 
M1 + M2 

M1M2 
M1 + M2' 



r rO 



re 



w . f M1M2 

= wrliAMi + AM2) 
+ r2v{AMi + AM2) 
M1M2 



re 



MI + M2 



(12) 



thus 
A{re) = 
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The Ah/h term is: 
Ah _ AMi AM2 
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2-he2 A(Mi-FM2) V 



2(1 -e2)5 



Mo 
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Ml (Ml + M2) ' M1+M2 

Combining Eqs. (9), (10), (15) and (17), we can obtain 
the changes of orbital semi- major axis A and eccentricity 
e: 
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Fig. 4. The fitting of the predicted to observed heavy- 
element abundances of 14 barium stars in standard case 
of wind accretion. But the curves in Fig. 4j and k represent 
the results of higher accretion rate: 0.5 times of the Bondi- 
Hoylc's rate. In every figure, the alphabet 'a' represents 
the times of the corresponding neutron exposure in the 
^^C profile suggested by the Fig. 1 of GaUino et al. (1998), 
and 'P' represents the orbital period of the barium star. 



